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Progress to Date 


This report documents the progress that has been made in the proposed research work since the 

first annual report. The basic goals for year two have been completed. Specific 

accomplishments for the past year are itemized below. 

1. A paper has been presented that documents the work using the T-501 engine model. The 
paper was presented at the Symposium on Nonlinear and Stochastic Dynamics, held at 
the 1994 ASME Winter Annual Meeting, November 13-18, 1994, Chicago, Illinois. The 
paper has also been submitted to the ASME Journal of Vibration and Acoustics. A copy 
of the paper is included in the Appendix. 

2. A paper that discusses some of the work using the T-501 engine model has been 
presented at the First Industry /University Symposium on High Speed Civil Transport 
Vehicles, December 4-6, 1994, Greensboro, North Carolina. A copy of the paper is 
included in the Appendix. 

3. Development of the experimental facilities has continued. This includes a new design 
for the magnetic bearing and the inclusion of housing dynamics in an additional test rig. 

4. A paper that documents an experimental/simulation study of auxiliary bearing 
rotordynamics has been presented at the 40 h ASME International Gas Turbine and 
Aeroengine Conference, June 5-8, 1995, Houston, Texas. A copy of the paper is 
included in the Appendix. 

5. A paper that describes a rotordynamical model for a magnetic bearing supported rotor 
system, including auxiliary bearing effects has been written and submitted to the ASME 
15th Biennial Conference on Vibration and Sound, to be held in Boston, Massachusetts, 
Sept. 17-21, 1995. A copy of the paper is included in the Appendix. 

6. A finite element model for a foil bearing has been developed. Rotordynamical studies 
of a rotor supported by foil bearings is currently underway. This includes a study of the 
influence of nonlinear effects arising from foil bearings on rotordynamic behavior. 

7. Additional studies of rotor/bearing/housing dynamics are currently being performed with 
the experimental rig and simulation models. 

8. The effects of sideloading on auxiliary bearing rotordynamics is being studied using the 
magnetic bearing supported rotor model. 

9. Two students affiliated with this project have graduated with M.S. degrees - April M. 
Free and James L. Lawen, Jr. 
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ABSTRACT 

This paper investigates the steady-state responses of a ro- 
tor system supported by auxiliary bearings in which there 
is a clearance between the rotor and the inner race of the 
bearing. A simulation model based upon the rotor of a pro- 
duction jet engine is developed and its steady-state behavior 
is explored over a wide range of operating conditions for var- 
ious parametric configurations. Specifically, the influence of 
rotor imbalance, support stiffness and damping is studied. 
It is found that imbalance may change the rotor responses 
dramatically in terms of frequency contents at certain oper- 
ating speeds. Subharmonic responses of 2nd order through 
10th order are all observed except the 9th order. Chaotic 
phenomenon is also observed. Jump phenomena (or double- 
valued responses) of both hard-spring type and soft-spring 
type are shown to occur at low operating speeds for sys- 
tems with low auxiliary bearing damping or large clearance 
even with relatively small imbalance. The effect of friction 
between the shaft and the inner race of the bearing is also 
discussed . 

NOMENCLATURE 

Cb — auxiliary bearing support damping, lb.s“/in. 

Cb\1> — auxiliary bearing torsional damping, Ib.in.s 
F n = normal force, lb 
F t = friction force, lb 

F x = external force vector acting on the rotor in X direction 
Fy = external force vector acting on the rotor in Y direction 
I a = rotor inertia matrix 


J b = moment of inertia of auxiliary bearing, Ib.in.s 2 
Kb = auxiliary bearing support stiffness, lb/in. 

Kc = contact stiffness, Ib/in. 

Mb — auxiliary bearing mass, lb.s 2 /in. 

Mk = mass of kth rotor element, lb.s 2 /in. 

N = total number of modes considered 
N B 1 = node number at auxiliary bearing #1 
N B2 = node number at auxiliary bearing # 2 
Qx = rotor modal coordinate vector in X direction 
Qy = rotor modal coordinate vector in Y direction 
Rq = radius of auxiliary bearing bore, in. 

Rm = radius of auxiliary bearing pitch, in. 

Rb = radius of rotor journal, in. 

Xh =s rotor physical coordinate vector in X direction 
Y/* = rotor physical coordinate vector in Y direction 
e = rotor imbalance eccentricity, in. 
g = gravitational acceleration, in./s 2 
t = time, s 

A = deformation at the contact point, in. 

r = 

= rotor free-free modal rotation matrix 
Q = rotor operating speed, rad/s 
$ = rotor free-free modal displacement matrix 
8 = Rb — Rr , auxiliary bearing clearance, in. 

= dynamic friction coefficient 
ji ^ = roiling friction coefficient 

1 1>Q = angular displacement of auxiliary bearing inner-race 
£ = modal damping coefficient 
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INTRODUCTION 

One of the most innovative developments in the turbo- 
machinery field involves the use of active magnetic bearings 
(AMB) for rotor support. This technology provides the po- 
tential for significant improvements in the dynamic behavior 
of rotor systems, allowing for loading, eccentricity, shaft po- 
sition and vibration to be continuously monitored and con- 
trolled. In order to protect the soft iron cores of the magnetic 
bearings and to provide rotor support in the event of failure 
of the bearing or during an overload situation, backup (or 
auxiliary) bearings, with a clearance between the rotor and 
the inner race of the bearing, are usually included in the ro- 
tor design. This clearance introduces a nonlinear dynamical 
feature which may significantly impact the behavior of the 
rotor. 

Magnetic bearing systems appear to provide particularly 
great promise for use in aerospace applications. There are 
active programs at many of the major jet engine manufac- 
turers to develop engines supported by magnetic bearings. 
Safety is a major concerns in any aeronautical design. To- 
ward this end, it is desirable to design the rotor system to 
take maximum advantage of the backup bearings and use 
them as true auxiliary bearings to provide support during 
critical situations in a safe and consistent manner. An im- 
portant concern in this regard is the dynamic behavior of the 
rotor when it comes into contact with the auxiliary bearing. 
If safe and effective operation of the engine is to be ensured 
during these periods, it is essential that designers have a very 
good understanding of the steady-state dynamics of rotor 
systems with clearance effects. 

There are a number of studies in the literature concerned 
with the dynamics of rotors with clearance effects. Ya- 
mamoto (1954) conducted a systematic study of rotor re- 
sponses involving bearing clearance effects. Black (1968) 
studied the rotor/stator interaction with a clearance. He 
concluded that rotor/stator interactions may occur in a vari- 
ety of forms and circumstances, including jump phenomena. 
Ehrich (1966) reported the first identification of a second 
order subharmonic vibration phenomenon in a rotor system 
associated with bearing clearance (1966). Bently (1974) pub- 
lished experimental observations of second and third order 
subharmonic vibration in a rotor system. Later, Muszynska 
(1984) cited the occurrence of second, third, and fourth order 
subharmonic responses in a rotor rubbing case and Ehrich 
(1988 and 1991) observed eighth and ninth order subhar- 
monic vibration as well as chaotic vibration in a high speed 
turbomachine. Childs (1979 and 1982) published two papers 
to explain the mechanism for the second and third order 
subharmonic responses noted above. He stated, with great 
insight, that a motion due to nonsymmetric clearance effects 
is a fractional-frequency phenomenon.” 

While those studies have greatly enhanced the understand- 



FIG. 1 DIAGRAM OF THE FEM ROTOR MODEL 


ing of clearance effects on rotor dynamics, a more detailed 
understanding of the dynamical behavior of such systems is 
needed. The perspective of much of this earlier work is that 
the clearance exists as a result of manufacturing error or mis- 
fitting. That is, it is due to an abnormal situation. However, 
in a rotor system fitted with magnetic bearings and auxiliary 
bearings, the clearance becomes a design parameter rather 
than an irregularity. From this point of view, it is impor- 
tant to develop a detailed quantitative understanding of the 
dynamic responses that are to be expected. Such knowledge 
will provide guidelines for the selection of auxiliary bearing 
parameters. 

It seems that there have been little work to date that is 
specifically concerned with auxiliary bearings in magnetic 
bearing supported rotor systems. Two papers that are di- 
rectly related to research on auxiliary bearings were both fo- 
cused on transient responses. Gelin et ai. (1990) studied the 
transient dynamic behavior of rotors on auxiliary bearings 
during the coast down. Ishii and Kirk (1991) investigated 
the transient responses of a flexible rotor during the rotor 
drop after the magnetic bearings become inactive. In both 
papers, idealized rotor models are used and it is assumed 
that once the magnetic bearings fail, the torque is cut off 
and consequently the rotor speed approaches zero. 

In this paper, simulation results are presented for a 
complex rotor system supported by auxiliary bearings with 
clearance at each end of the rotor. This work is specifically 
concerned with systems in which the clearances are quite 
small (on the order of a few mils), which is appropriate 
for jet engine applications in which the backup bearing is 
acting to provide rotor support on a consistent basis. The 
influence of rotor imbalance, support stiffness and support 
damping are investigated using direct numerical integration 
of the governing equations of motion and the harmonic 
balance method. Some insights are obtained with regard 
to the frequency and amplitude behavior of the steady-state 
vibration of such a system. 
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coordinate transformation: 



FIG. 2 AUXILIARY BEARING MODEL 

SIMULATION MODEL 

The rotor is modelled using free-free normal mode shapes 
and natural frequencies obtained through finite element 
analysis. The model data is representative for the rotor 
of a jet engine. Fig. 1 shows a schematic diagram of the 
FEM rotor model. Parametric information about the model 
is listed in Table 1. The torsional motion of the shaft is not 
considered in this paper. Using state space representation 
and modal coordinates, the equations of motion for the rotor 
are expressed as 


jV 

X Bk = ^22 $1* Qx%> 
1=1 
N 

Yuk = ^2 ® ik Q yi ' 
»= i 


(k = NB1,NB2) 


The equations of motion for the auxiliary bearings are 
derived using the model shown in Fig. 2 


M B kX B k + CskXsk 4 - KskXsk 

= F nk cos a* - Ftk sin a k + M Bk g, (2. a) 
MskYBk + C Bk Y B k 4- FBkYBk 

= F nk sin a* + F tk cosat, (2.6) 

JBkVBk + C B yi>k^Bk 

= F t kR B k — HrpFnkRmk* (2 .c) 


where 


a* = tan 1 


Y/u - YBk 
Xftk - X B k 


(k = NBl y NB2) 


At this point, the rotor and the back-up bearings appear 
to be uncoupled. However, the force vectors Fx and Fy 
on the right hand sides of equations (1) are partially due to 
rotor/auxiliary bearing interaction. In fact, we have 


Q*+ 2 Cu; n Qx 4 . QTQy + 

+ 2fXw n Qy = & t Fx, 
Qy-b 2 Cw„Qy - QTQx + 
-2QC^nQx =^ r Fy, 


where 

Fx = Fxm }~ 1 y 

Fy = {Fy 1 , Fy 2> . . . , Fy m }~ l y 

Qx =#-%, 

Qy 


Fxk = - F nk cos a k + Ftk sin a k 
+ M k g + M k e Q 2 cos 
Fy k = - F nk sina k - F tk cosa k 
+ M k e ft 2 sin (Qt). 


The rotor/bearing interaction is represented with the normal 
force F nk 


where 


r K c 6 k , A k < 0, 
1 0 , A*>0, 


(3.a) 


with 

x« = {Xfti, Xr2, . . . , vYflm} -1 

(m = total number of nodes) 


Ak = (XRk — X B k)cos Q k 4* (V/u - YBk) sin a k - 6 k 

and the Coulomb friction force Ftk- As long as there exists 
slip at the contact point, the friction force obeys 

F lk = fkF nk . (3.6) 


The physical displacements of the rotor at the two aux- 
iliary bearing locations can be obtained using the following 
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However, when there is no slip at the contact point, the 
friction forces are solved from equations (1) and (2) using 



the kinematic constraint that the circumferential velocities 
of the rotor and the inner-race of the back-up bearing at 
the contact point equal to each other. At the same time, if 
this solved friction force exceeds the maximum static friction 
force (= li s F n k), equation (3.b) applys again. 

DISCUSSION OF RESULTS 

The rotor is modeled with 34 stations (as shown in Fig. 1) 
and the first four modes (two rigid body and two flexible 
modes) are included in the simulation model. The two 
auxiliary bearings are located at nodes 3 and 33, respectively. 
This arrangement is taken to represent one of the most 
technically feasible configurations in that it greatly simplifies 
bearing maintenance. It is assumed that the two auxiliary 
bearings are identical in terms of stiffness, damping and 
friction characteristics. Some nominal system parameters 
used for the simulation study are £=0.03, 

R m j c = l.l Rbic, /i 5 = 0.5 and ^=0.002. To avoid excessive 
cluttering of plots, all the results that are presented in this 
paper correspond to node 3, the location of bearing #1. 

Since the total system which includes two bearings and 
associated friction forces as well as the inner-race motions 
is rather complicated and requires considerable amount of 
computer time for the solutions to converge, the friction 
effect is examined first to see if the model can be further 
simplified. It turns out that the steady-state results obtained 
with and without friction are virtually identical. Even the 
differences in transient responses are quite small, as can be 
seen in Fig. 3(a) and 3(b). The only remarkable effect is 
on the transient responses of the inner-race as shown in 
Fig. 3(c). This observation is confirmed by numerous runs 
using different system parameters and rotor speeds. The 
lack of significance of friction may be attributed to several 
factors. First, the inertia of the inner-race is quite small in 
comparison to the rotor mass, so vibration of the bearing 
has little influence on the rotor vibration. Second, ball 
bearings exhibit quite negligible torsional resistance under 
normal conditions. As a result, the terms that are related 
to the friction forces and the rotational motion of the inner 
race are not included in the simulations that are discussed in 
the following paragraphs. 

The steady-state response characteristics of the system 
are obtained through numerical integration of the simplified 
version of governing equations (1) and (2). Near-zero initial 
conditions are used, simulating situations where the AMBs 
are functioning properly prior to a system failure. Multiple 
solutions with other initial conditions are not sought at 
present. 

It is well known from linear analyses that imbalance 
greatly affects the steady state vibration amplitudes of a 
rotor system. However, it is observed from the current 
work that imbalance may also influence frequency content of 


the rotor responses quite dramatically at certain operating 
speeds. A typical case with such an imbalance effect is shown 
in Figs. 4 and 5, where orbits and corresponding frequency 
spectra of the rotor for different values of imbalance at the 
speed of ft = 1000 are plotted. For this particular case, 
there exist eight ranges of imbalance values that result in 
eight different types of rotor responses. 

For e < 0.0009, the rotor rotates near the bottom of 
auxiliary bearings and the responses are predominantly syn- 
chronous (Figs. 4(a)-4(b) and 5(a)-5(b)). As imbalance in- 
creases, the 2ft superharmonic component approaches the 
magnitude order of the synchronous component (Figs. 4(b) 
and 5(b)). However, the responses are of small amplitude. 
For 0.0010 < e < 0.0013, the responses are dominated by 
ft/2 subharmonic components (Figs. 4(c) and 5(c)). In other 
words, the amplitude of fl/2 component is greater than that 
of the synchronous. For 0.0014 < e < 0.0016, the ft/2 
subharmonics disappear and the ft/3 subharmonics become 
dominant (Figs. 4(d) and 5(d)). So far, the overall amplitude 
of the responses are not large, the rotor just bounces near the 
bottom of the auxiliary bearings. For 0.0017 < e < 0.0027, 
the orbits become chaotic-looking (Fig. 4(e)) and the spec- 
trum contains a lot of noise (Fig. 5(e)). In this range of 
imbalance, the rotor changes from bouncing near the bot- 
tom to bouncing around the full clearance of the bearing as 
imbalance increases. In the middle of this transition range, 
true chaos is observed. The Poincare map shown in Fig. 6(a) 
and the frequency spectrum shown in Fig. 6(b) demonstrate 
that the response has all the characteristics of a chaotic 
phenomenon. It should be noted that even though the or- 
bits are chaotic looking, the amplitudes are not the largest 
among all the cases for this particular parametric configu- 
ration. For 0.0028 < e < 0.0034, the orbits orbits are no 
longer chaotic-looking (Fig. 4(f)). The spectrum shows they 
are ft/5 subharmonic responses (Fig. 5(f)). Notice the am- 
plitudes are the largest for this parametric configuration. For 
e = 0.0035, the amplitude suddenly becomes smaller even 
though the imbalance has become larger (Fig. 4(g)). And 
the rotor bounces near the bottom of the auxiliary bearings 
again. The frequency spectrum shows it is ft/8 subharmonic 
response (Fig. 5(g)). For 0.0036 < £ < 0.0042, the orbits 
become chaotic-looking again (Fig. 4(h)). But the frequency 
spectra are very similar to the subharmonic cases (Fig. 5(h)), 
only with some discrete noise. Finally, for e > 0.0043, the re- 
sponses become predominantly synchronous again (Figs. 4(i)- 
4(j) and 5(i)-5(j)). But this time as imbalance increases, the 
2ft superharmonic component become smaller and smaller 
(Figs. 4(j) and 5(j)). 

Examining all the orbits in terms of amplitudes as imbal- 
ance increases, we can see the characteristics of a jump-type 
phenomenon (Cunningham, 1958). The jump— down takes 
place around 0.0034 < e < 0035 where the rotor jumps 
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from full-clearance bouncing to near-bottom bouncing. Fur- 
ther investigation is needed to better understand this type 
of change- 

imbalance responses at some other operating speeds and 
for other parametric configurations exhibit similar changes 
as imbalance varies, though the imbalance ranges and cor- 
responding response types may not be as well defined as 
in the above cases. In fact, subharmonic responses from 
ft/2 through ft/ 10 are all observed except ft/9 as shown 
in Figs. 7. Surprisingly, those subharmonics are not di- 
rectly related to the system’s natural frequencies as were the 
cases with other researchers’ findings (such as Ehrich, 1988). 
Moreover, several types of subharmonic responses may occur 
at a single operating speed. It should be noted that Chen et 
al. (1993) also reported occurrence of three stable subhar- 
monic responses at a single rotor speed in a SFD supported 
rotor system but did not provide any explanation for their 
findings. In their case, even the imbalance did not vary. Ap- 
parently, further research is needed to find a mechanism to 
explain this multi-subharmonic vibration phenomenon. On 
the other hand, it should be pointed out that these subhar- 
monic responses are not typical cases. While some of them 
are observed to exist within a certain range of parameters, 
the majority of them occur only for some specific parametric 
configurations. 

Due to space limitations, results for other parametric 
configurations are not systematically plotted. A general 
summation of the observations is presented instead. A 
common feature among all the responses is that for very 
small imbalance, the responses are always synchronous. 
The imbalance range that result in synchronous dominated 
responses depend on several system parameters. For small 
back-up bearing stiffness (such a s Kq — Q.213e+6) and 
normal damping ( Cq = 157.0), the responses are almost 
always synchronous. Only ft/2 subharmonic are observed 
at a few operating speeds with a very narrow range of 
imbalance. It should be noted that even though a lower 
Kb may leads to a better system response, it may also 
fail to protect the magnetic bearings due to the fact that 
it could result in a larger rotor orbit-center offset. However, 
the dramatic response changes discussed above may occur 
again if the damping becomes small (such as Cq = 57.0) 
even though the stiffness still remains small. On the other 
hand, increasing the damping Cq alone may not be able 
to eliminate those dramatic changes. It is observed that 
those changes can still occur for Cb being as large as 
700.0. Reducing the size of clearance S may not eliminate 
the response changes at certain speeds. But it can narrow 
the operating speed range where those changes occur. For 
example, response changes are eliminated for ft > 1500 
when 6 is reduced from 0.002 to 0.001 with all other 
parameters remaining the same, but response changes still 


occur for ft < 1400. 

It is obvious that rotor responses involving nonsymmetric 
bearing clearance effect are very complex problems and 
numerical integration alone is not a sufficient tool to obtain 
a global picture of the system responses. The harmonic 
balance method is then used for the investigation of global 
system behavior. However, it is only attempted for situations 
with very small imbalance values. In addition, only the lft 
harmonic is considered. The complex frequency contents 
associated with medium and large imbalance values makes 
it a formidable task to apply the harmonic balance method 
for other cases. Nevertheless, some useful information can be 
drawn from these results. After all, an adequately balanced 
rotor system should have very small imbalance under normal 
conditions. 

Fig. 8(a) shows that nonsymmetric clearance effect is 
equivalent to asymmetric stiffness effect with regards to 
critical speeds. The clearance actually splits the first critical 
speed into two pseudo-critical speeds. In the X direction, 
the gravity force tends to keep the rotor in contact with 
the bearing at low operating speed. Thus, the apparent 
stiffness is almost the same as Kb a>nd the pseudo-critical 
speed is nearly the same as the critical speed for the linear 
case (<5 = 0). But in the Y direction, the clearance 
results in a lower apparent stiffness and, consequently, an 
additional lower-value pseudo-critical speed. It is seen that 
several higher order additional pseudo-critical speeds are 
created in the operating speed range in addition to the 
1st additional pseudo-critical speed. It is noted that the 
response in the X direction also departs from the linear 
case at high operating speed. This is because the imbalance 
force becomes dominant at high rotor speed which in turn 
makes the gravity force less significant and the clearance 
effect more important. Fig. 8(b) shows that changing the 
auxiliary bearing stiffness has little effect on the pseudo- 
critical speeds of the system. However, for a larger value 
of imbalance, a higher Kb does leads to a greater tendency 
of double— valued responses. In each direction, for either a 
stiffness increase or an imbalance increase, the 1st pseudo- 
critical peak tend to become a hard-spring type jump and 
the 2nd one tend to develop into a soft-spring type jump, 
with the tendency decreasing as the pseudo-critical’s order 
increases. 

Fig. 9(a) shows double-valued responses in the Y direction 
for four different values of clearance. It is seen that a larger 
clearance results in wider rotor speed range of double-valued 
responses. It is also observed that as clearance increases, 
the apparent stiffness decreases and the first pseudo— critical 
speed shifts to a lower value. Fig. 9(b) shows the double- 
valued responses in the X direction. Even though the jumps 
themselves are smaller in magnitude, they are more obvious 
in trend. Notice how little the change is for the first pseudo- 
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critical speed in the X direction. 

Fig. 10(a) shows the influence of auxiliary bearing damping 
on the double-valued responses in the Y direction. It is 
observed that the damping has to be quite large to eliminate 
the double-valued responses associated with the first pseudo- 
critical speed. Fig. 10(b) shows the influence of Cb on tKe 
double-valued responses in the X direction. In both figures, 
it should be noted that as Cb decreases, the second pseudo- 
critical speed peak will develop into a soft^spring type jump 
and the third pseudo-critical speed peak will evolve into a 
hard-spring type jump. 

The system behaviors for higher operating speed range are 
not shown in Figs. (9) and (10) so that the jump phenomena 
can be more clearly illustrated. It is also because that the 
system’s responses at high operating speed range with the 
same parameters are more or less regular, in other words, 
mainly amplitude changes. 

CONCLUSIONS 

As a summary of the results discussed above, the following 
conclusions can be drawn: 

1. Imbalance may change the rotor responses dramatically in 
terms of frequency contents at certain operating speeds, 
especially under conditions of large clearance, high bearing 
stiffness and low bearing damping. With imbalance 
changing, as many as eight different types of responses 
may occur for a particular parametric configuration at a 
single operating speed. 

2. Subharmonic responses of second order through tenth 
order are all observed except for the ninth order case. 
However, the majority of them are not typical cases, 
and were observed only for quite particular parametric 
configurations. 

3. Chaotic phenomenon is observed to occur occasionally. 
However, the amplitudes associated with such motion are 
not among the largest. 

4. Nonsymmetric clearance effects influence the critical speeds 
in a manner similar to asymmetric stiffness effects. 

5. Double-valued responses in the form of both hard-spring 
type jump and soft-spring type jump are observed to be 
possible at low operating speeds with low auxiliary bearing 
damping or high imbalance. With large clearances or high 
bearing stiffness, the jump phenomenon may occur for 
even relatively small imbalances. 

6. The effect of friction between the shaft and the inner race 
of a rolling element auxiliary bearing on the dynamics of 
the rotor is quite small and can reasonably be neglected 
for steady-state analyses. 
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TABLE T FEM MODEL DATA 


Station 

O.D. 

I.D. 

\VT. 

Polar Mom. 

(in.) 

(in.) 

(in.) 

(ib) 

(Ib.in.s 2 ) 

5.400 
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2.000 

0.000 

0.000 
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2.000 
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2.000 
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0.000 

6.505 

7.300 

6.950 

0.100 

0.100 

7.100 

7.300 
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8.670 

7.300 

7.120 

30.400 

425.000 

10.900 

8.400 

8.220 

0.000 
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13.120 

9.300 

9.120 

0.000 

0.000 

15.310 

10.100 

9.920 

20.300 

538.000 

17.280 
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9.150 
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0.000 
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12.150 
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23.670 

12.400 
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24.900 

12.500 

12.350 
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26.120 
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12.350 
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27.400 

12.500 

12.350 

22.550 ! 

539.000 

27.750 

12.500 

2.500 

1.027 
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2.500 
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25.022 

49.761 

37.000 
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17.961 

0.000 
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4.500 
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17.961 

0.000 

18.000 
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2.500 

16.328 

0.000 
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2.500 

16.328 

0.000 

58.000 

4.500 

2.500 

14.614 

0.000 

61.950 

4.500 

2.500 

6.531 

0.000 

62.000 

6.000 

2.500 

0.954 

0.000 

62.405 

6.000 

5.670 

34.935 

1216.500 

64.600 

6.000 

5.670 

0.000 

0.000 

67.000 

6.000 

5.670 

117.000 

3750.500 

(59.500 

6.000 

5.670 

0.000 

0.000 

(59.505 

2.850 

1.250 

0.000 

0.000 

72.500 

2.850 

1.250 

10.197 

12.000 

73.600 

2.850 

1.250 

0.000 

0.000 


Modulus of elasticity: E=2.80e-f 7 psi 
Shear Modulus of elasticity: ( l .OSe-t-7 psi 
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FIG 3 EFFECTS OF BEARING FRICTION (C s = 157.0, 
/v'a=0.313e+6. 0=1400. 6=0.002. e=0.0002) 
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FIG. 4 IMBALANCE RESPONSES - ORBITS 
(<5=0.002, A'*=0.313e+6 t C fl = 157.0) 
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FIG. 7 SUBHARMONIC RESPONSES (ORBITS) 
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FIG. 8 INFLUENCE ON CRITICAL SPEEDS FIG. 10 EFFECTS OF BACK-UP BEARING DAMPING 

(a) e=0.0001, Cz? = 157.0, 7v B =0.313e+6 (<5=0.002, e=0.0001, A' s =0.313e+ 6) 
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FIG. 9 CLEARANCE EFFECT ON JUMP PHENOMENA 
(e=0.0001, A'5=0.313e-h6, Cij=300.0) 


FIG. 6 CHAOTIC RESPONSES (e = 0.0023, 
6=0.002. A'5=0.313e+6, C fl = 157.0) 
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ABSTRACT 

This paper presents a study of the dynamic behavior of a rotor system supported by 
auxiliary bearings. The steady - state behavior of a simulation model based upon a 
production jet engine is explored over a wide range of operating conditions for varying 
rotor imbalance, support stiffness and damping. Interesting dynamical phenomena, 
such as chaos, subharmonic responses, and double-valued responses, are presented and 
discussed. 

KEYWORDS: magnetic bearings, auxiliary bearings, nonlinear rotordynamics. 

NOMENCLATURE 

Cb = auxiliary bearing support damping, lb. s/in. 

F n = normal force, lb 

Fx* = external force vector acting on the rotor in X direction 
Yy = external force vector acting on the rotor in Y direction 
I a = rotor inertia matrix 

I(b — auxiliary bearing support stiffness, Ib/in. 

Kc = contact stiffness, lb/in. 

Mb = auxiliary bearing mass, lb.s“/in. 

Mk = mass of kth rotor element, lb.s~ /in. 

N = total number of modes considered 
Nbi = node number at auxiliary bearing #1 
N B2 = node number at auxiliary bearing #2 
Qx = rotor modal coordinate vector in X direction 
Q y = rotor modal coordinate vector in Y direction 
R b = radius of auxiliary bearing bore, in. 

Rft = radius of rotor journal, in. 

Xk = rotor physical coordinate vector in X direction 
Yfi = rotor physical coordinate vector in Y direction 
e = rotor imbalance eccentricity, in. 
g = gravitational acceleration, in./s** 
t — time, s 

A = deformation at the contact point, in. 

r = ^ T i a ^ 

* = rotor free-free modal rotation matrix 
fi = rotor operating speed, rad/s 
<|> = rotor free-free modal displacement matrix 
6 = Rb - Rr, auxiliary bearing clearance, in. 


aECEDHMG PAGE BLANK NOT FiLMEb 
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( = modal damping coefficient 


INTRODUCTION 

Active magnetic bearings are one of the most innovative recent developments in the turbomachincry 
field. This technology provides the potential for significant improvements over other types of rotor 
support, including elimination of wear and bearing friction-related energy losses as well as a means 
of actively suppressing rotor vibration. A critical component of any magnetic bearing design is the 
auxiliary bearing, which protects the soft iron core of the magnetic bearing and provides rotor support 
in case of overload or failure of the primary (magnetic) bearing. Magnetic bearing systems appear 
to provide particularly great promise for use in aeronautical applications. In this regard, current 
effort is directed toward developing jet engines with rotors supported by magnetic bearings. For such 
applications, safety is an important concern. Toward this end, it is desirable to design the system to 
operate with auxiliary bearing support for an extended period of time. 

A number of different bearing types have been suggested as auxiliary bearings. These include bushings, 
rolling element bearings, and various types of journal bearings. The most commonly considered are 
rolling element bearings. The major disadvantage associated with using rolling element bearings (or 
bushings) is the requirement of a clearance between the rotor and the inner race of the bearing, without 
which many of the advantages associated with using magnetic bearings would be reduced or eliminated. 
This clearance introduces a nonlinear dynamical feature which may significantly impact the behavior 
of the rotor. 

There are quite a number of studies in the literature concerned with nonlinear rotordynamics. 
Yamamoto (1954) conducted a systematic study of rotor responses involving bearing clearance effects. 
Black (1968) studied the rotor/stator interaction with a clearance. Ehrich (1966, 1988 and 1991), Bently 
(1974), Muszynska (1984) and Childs (1979 and 1982) observed and studied subharmonic responses 
associated with clearance effects. While this work has served to greatly enhance the understanding 
of such systems, more detailed study is needed. Much of this earlier work was conducted from the 
perspective that the clearance is a result of manufacturing error or misfitting and is best eliminated. 
However, in a rotor system fitted with magnetic bearings and auxiliary bearings, the clearance becomes 
a design parameter rather than an irregularity. From this point of view, it is important to develop a 
better understanding of the expected dynamic responses. Such knowledge will provide guidelines for 
the selection of auxiliary bearing parameters. 

There has been relatively little work available in the open literature that is specifically concerned with 
auxiliary bearings. Two papers that are directly related to research on auxiliary bearings are Gelin 
et al., (1990) and Ishii and Kirk (1991). Both of these papers are mainly concerned with transient 
responses. The current work is specifically concerned with developing a better understanding of the 
expected steady-state dynamical behavior for an auxiliary bearing supported rotor. Simulation results 
for a rotor based upon a production jet engine are presented and discussed. 

The model used in the current study has two principal components - the rotor and the auxiliary 
bearings. The rotor is modelled using the free-free bending mode shapes and natural frequencies 
obtained through finite element analysis. The finite element code uses 31 stations and the first four 
modes (two rigid bodv and two flexible modes) are included in the simulation model. Figure 1 shows 
a schematic diagram of the finite element model. The rotor equations of motion can be expressed in 
terms of modal coord i nates as: 


Q.Y+2Cw n Q.v + erQr+w'Q.Y 

+ 2nC w nQl r = ^ r F.Y, 0- a ) 

Qr+2<WnQ>- - nrQ.Y - --'-Qy 
- 2Q(u n Q x = <& r F r . 


( 1 . 6 ) 
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FIGURE 1 DIAGRAM OF THE FINITE FIGURE 2 AUXILIARY BEARING MODEL 

ELEMENT MODEL 


where 


with 


F.V = {F.V1, F.Y2, • 

• • , Fxm}~ T 

Ol 

£ 

£ 

11 

fa 

■ , FYm}~ T ■ 

Qx = ^ l X fl , 


Q Y = ^"‘Yr, 



/Wn, 0 

... 0 x 

u; n = 

0 

... 0 


\ 0 0 

... U> nN ) 

Xft = {A'ri, A'fl2, . 

...x Rm }- T 

Yr = {Vm , Vri, . . 

"Y Rm }~ T . 


(m = total number of nodes) 


The physical displacements of the rotor at the two auxiliary bearing locations can be obtained using 
the following coordinate transformation: 


.V 

X Rk = ^ $ki Qxi, 

N 

Y Rk = ^2 Qyu 

i=t 


(k = i'V 01> iV 02 ) 


The model for the auxiliary bearings is shown in Figure 2. The governing equations of motion are: 


Bk + C Bk A Bk + A SfcA 

= Fn^cosa* + jV/ 0fc g, (2.a) 

MskYBk + C Bk Yq k + A Bk^ Bk 

= ^ aina*, 


( 2 . 6 ) 
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where 


Qjt = tan 1 


Yttk ~ }j3k 

X Rk - A'flfe ‘ 


(& = A bi, X 32) 


Fxk = - F n k COS a k + M k g + M k e ft 2 cos {Sit), 
Fy k ~ F nk sin a k + M k e fi 2 sin (Qt). 

The rotor/bearing interaction is represented with the normal force F nk 

j? _ j K c S k , A k > 0, 
n * - \0, A k < 0, 

where 

Afc = (Xflfc - A"bA:) C05 a/; + (Yflfc - YBk) sin a h - $lc 


( 3 ) 


DISCUSSION OF RESULTS 

The model described in the previous section was used to study the steady-state dynamic behavior 
of such a system using direct numerical integration of the equations and a harmonic balance code. 
Studies were performed for varied values of imbalance, support stiffness, and support damping, 
respectively. For the purposes of the current work, the two auxiliary bearings are located at nodes 
3 and 33, respectively. It is assumed that they are identical in terms of stiffness, damping and 
friction characteristics. The nominal system parameters used for the simulation study are Mbi =0.0023, 
Mb2 = 0.0024, A r c=‘2.855e+6, and £=0.03. All the results that are presented correspond to node 3, 
the location of bearing 1. 

The dynamical behavior of a rotor supported by bearings with clearance coupled with the nonsymmetry 
resulting from gravitational effects can be quite complex. The harmonic balance method is first used 
to investigate the synchronous behavior for small imbalance values. (Please note that the complex 
frequency contents associated with medium and large imbalance values makes it a formidable task 
to apply the harmonic balance method for other cases.) Figure 3 is a typical plot of the steady 
state response amplitudes as functions of the rotor speed with and without bearing clearance. The 
nonsymmetric effects resulting from gravity loading is similar to what occurs for a rotor supported by 
nonsymmetric bearing stiffnesses with regard to mlluence on critical speeds. The first critical speed 
splits into two distinct values. For the X direction, the gravity force tends to keep the rotor in contact 
with the bearing at low operating speed. The apparent stiffness in this direction is approximately the 
same as Kq and the resulting critical speed about the same as the critical speed for the linear case 
(£=()). For the Y direction, the clearance results in a lower apparent stiffness and a lower critical speed. 
Several higher order pseudo-critical speeds are also created in the operating speed range (about 1500 
rad. /see). It should be noted that t lie response in the X direction also departs from the linear case 
at high operating speed. This is because the imbalance force becomes dominant at high rotor speed 
which in turn makes the gravity force less significant and the clearance effect more important. 

Figure 4(a) shows some typical results for varying imbalance. Some of the more dynamically interesting 
results occur for cases of relatively large imbalance. It is observed that imbalance mav influence the 
frequency content of the rotor responses quite dramatically at certain operating speeds. There exists as 
many as eight ranges of imbalance values that result in eight different types of rotor responses. In fact, 
subhariiionic responses from f)/‘2 through 11/10 are observed. Those subharmonics are not directly 
related to the s\ stem’s natural frequencies as were t he cases with other researchers' findings (such as 
Ehrich, 1988). Moreover, several types of subharmonic response.- may occur for identical parametric 
configurations, but different imbalance values. 

Figures 4(b) and 4(c) show typical results using Kb or C u as the variable parameter. Clear routes 
to chaos are observed. A- Kb increases beyond certain value, a period-doubling bifurcation always 
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FIGURE 3 INFLUENCE OF CLEARANCE 
ON CRITICAL SPEEDS 
(e=0.0001 , C b = 150.0, 
K b =0.313x10 6 ) 


FIGURE 4(a) BIFURCATION DIAGRAM 

WITH IMBALANCE AS THE 
VARIABLE PARAMETER 
(11=1000, <5=0.002. Cb = 150. 
K b =0.313x10 6 ) 


takes place. As Cg decreases below certain value, the responses always become chaotic. However, the 
bifurcation type is not well defined as for varied Kb- It should be pointed out that even though a lower 
Kg may leads to a better system response, it may also fail to protect the magnetic bearings due to the 
fact that it could result in a larger rotor orbit-center offset. In addition, higher bearing damping does 
not necessarily result in synchronous responses. It is observed that, with small imbalance (e < 0.0005), 
the responses are always synchronous for all speeds in the speed range considered (< 1800) if Cg is 
large enough (Cs = 250.0) and the clearance is small enough (<5 < 0.001). As with the other cases, 
higher bearing damping or lower bearing stiffness or lower imbalance tends to increase the probability 
of synchronous responses. 
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0.0010 
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FIGURE 4(b) BIFURCATION DIAGRAM 

WITH BEARING STIFFNESS 

AS THE VARIABLE 

PARAMETER 

(11=1500, <5=0.002, C e = 150, 

e=0.0008) 


FIGURE 4(c) BIFURCATION DIAGRAM 
WITH BEARING DAMPING 
AS THE VARIABLE 
PARAMETER 
(11=1500, <5=0.002, e=0.0015, 
K b =0.313x10 6 ) 
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CONCLUSIONS 

As a summary of the results discussed above, the following conclusions can be drawn: 

1. Imbalance may serve to dramatically alter the frequency contents of the rotor responses at certain 
operating speeds. This is particularly evident for cases of large clearance, high bearing stiffness and 
low bearing damping. 

2. For sufficiently high imbalance: 

i. There is a threshold level of damping below which complex dynamical behavior can be expected. 

ii. There is a threshold level of stiffness above which complex dynamical behavior can be expected. 

3. Clear routes to chaos are observed. As bearing stiffness increases beyond a certain vaiue, a period- 
doubling bifurcation takes place which leads to chaos. As bearing damping decreases below a certain 
value, chaos also tends to occur but the bifurcation type is not so clearly defined. 
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ABSTRACT 

Auxiliary bearings are a critical feature of 
any magnetic bearing system. They protect the 
soft iron core of the magnetic bearing during 
an overload or failure. An auxiliary bearing 
typically consists of a rolling element bearing 
or bushing with a clearance gap between the 
rotor and the inner race of the support. The 
dynamics of such systems can be quite complex. 
It is desired to develop a rotor dynamic model 
which describes the dynamic behavior of a 
flexible rotor system with magnetic bearings 
including auxiliary bearings. The model is 
based upon an experimental test facility. Some 
simulation studies are presented to illustrate 
the behavior of the model. In particular, 
the effects of introducing sideloading from the 
magnetic bearing when one coil fails is studied. 
These results are presented and discussed. 

NOMENCLATURE 

C = damping, N-s/m 

Cl = clearance of auxiliary bearing, m 

D = nominal gap thickness, m 

Fjr = external force vector acting on the rotor 

in 

X direction, N 

Fy = external force vector acting on the rotor 
in 

Y direction, N 

G = gravitational acceleration, m/s 2 


I = rotor inertia matrix 

t = current, amp 

K = stiffness, N/m 

k = gain value 

L = equivalent circuit length 

M = mass, kg 

m = total number of nodes 

N = total number of modes considered 

Nb i = node number at rotor left end 

Nb2 = node number at AMB 

jVj 93 = node number at auxiliary bearing 

Nimb = node number at imbalance location 

Nt = number of turns of wire per coil 

Q % = rotor modal coordinate vector in X 

direction 

Qy = rotor modal coordinate vector in Y 
direction 

Rb = radius of rotor journal, m. 

Rs = radius of auxiliary bearing bore, m. 
t = time, s 
t f = fail time, s 
v = voltage, volt 

tv = relative velocity at auxiliary bearing 
rotor/stator contact point 
X/t = rotor physical coordinate vector in 
X direction 

Xs = stator physical coordinate vector in 
X direction 

Y/t = rotor physical coordinate vector in 

Y direction 

Ys = stator physical coordinate vector in 

Y direction 

a = acceleration, rad/sec 2 
S = dynamic clearance 



r = * T i* 

$ = rotor free-free modal displacement matrix 
^ = rotor free-free modal rotation matrix 
-0 = imbalance vector 
fi 0 = permeability of free space 
p r = permeability of silicon steel 
= rotor operating speed, rad/s 
t = time constant 
6 = angular position of the shaft 
£ = modal damping coefficient 

Subscripts 
a = first 

B = auxiliary bearing 
b = second 
C = contact 
g — gain constant 
h = horizontal 
/ = left bearing 
r = rate constant 
v = vertical 
x = x — direction 
V = v — direction 
bi = bias 

1 = top magnet 

2 = right magnet 

3 =: bottom magnet 

4 = left magnet 

Superscripts 

a = current amplifier 
c = controller 
/ = lead 

p — sensor amplifier 
INTRODUCTION 

In recent years, the use of active magnetic 
bearings (AMB) for turbomachinery support 
has been an area of interest in both academia 
and in industry. Magnetic bearings provide 
the potential for significant improvements over 
other types of rotor supports, including elimi- 
nation of wear, bearing friction-related energy 
losses and a means of actively suppressing ro- 
tor vibration. A critical feature of any mag- 
netic bearing supported rotor system is auxil- 
iary bearings to protect the soft iron core from 
rotor contact (and subsequent damage) during 
an overload or failure of the AMB. The present 
work is concerned with developing a suitable 
model and investigating the dynamic behavior 
of such a system. 


EXPERIMENTAL TEST FACILITY 

The model and control system development 
presented in this study is for an experimental 
test facility. The apparatus consists of a radial 
magnetic bearing that is supporting the right 
end of a rotor. The left end of the rotor is 
supported by a ball bearing suspended in a 
frame by four springs. A photograph and the 
corresponding schematic diagram are shown in 
Figure 1. 

The apparatus has two basic components: 
the rotor shaft and the magnetic bearing. The 
shaft is made of steel and is 0.0098 meters in 
diameter and 0.4572 meters in length. A steel 
disk of diameter 7.62 cm, thickness 2.54 cm, and 
mass 0.810 kg is placed at the midpoint of the 
bearing span. Threaded holes on the disk allow 
for imbalance to be added to the system. The 
rotor is driven by a variable speed motor with 
a controller. Shaft vibration is measured using 
eddy current proximity displacement sensors 
fixed to measure displacement in both the 
vertical and horizontal directions. 



Figure l.a Experimental Test Facility 

The magnetic bearing is based upon a design 
described in Humphris, et al (1988). The basic 
parameters for the bearing are shown in Table 
1. The bearing consists of four electromagnets 
equally spaced around a soft -iron core. A pho- 
tograph and schematic diagram of the bearing 
are shown in Figure 2. 

SIMULATION MODEL 

A simulation model was developed lor the ro- 
tor system described above. The model has 
three principal components the rotor, mag- 
netic bearing, and electronics. 1 he governing 
equations for each are shown below. 





Figure l.b Schematic Diagram of Experimental 
Test Facility 



Figure 2, a Magnetic Bearing Assembly 



Figure 2.b Schematic Diagram of Magnetic 
Bearing 


The rotor is modelled using the free-free 
bending mode shapes and natural frequencies 
obtained through unite element analysis. The 
finite element code uses 11 stations and the first 
four modes (two rigid body and two flexible 
modes) are included in the simulation model. 
The rotor equations of motion can be expressed 
in terms of modal coordinates as 

Qx + flrQ„ + u>*Q* + $ T F r = 0, (l.o) 


Qy - firQx + u£Q v + * T Fy = 0, (1.6) 

The auxiliary bearing consists of a bail bear- 
ing with a clearance between the rotor and the 
stator. The governing equations are 

Bx^B "(■ KbxXb = Fi,auit(l.c) 
MbYb+CbhYb + I^By^B = Fy iau *, (1 -d) 

where 

FijJinC-^fcl) = ^Nn(^lxQr "b Cf r Qr) 
Fy f j« n (iVfci) = ^A^n(-^"/yQy "b C/yQy) 
Fx,imb(Nimb) = Sl 2 tp 6 cos6 + aip 6 sin9, 

Fy , im b(N imb ) = fl 2 Asin6 - afacosO, 

F x t aux(,N 53) = &N i 3 (FNCosf 3 - F Rsinf}) 
F yAUX (N b3 ) = * Nt3 (F N sinl3 + Y rcos/3) 
Y gr av = * T MG 
F n = <f>K c (6 - D). 

If v r = 0, then Fr = required static friction 
force to satisfy this condition. If v r ^ 0 or 
if required static friction force > fi s Fw, then 
F R = fJLk^N- 

F x = F x,lin 4* 4 “ F x?am 6 

F y = F y ^ n 4 F Vt tm6 4 F y|£lux - Fy fam 6 4 F grav 


where 








implemented using analog circuitry. The con- 
troller is modelled as 


Sin(/3) v/( X K - Xs ) 2 + (Yk - Y5) 2 , 

(a , Xfl-Xs 

cos( ^ ) v'CXr-x^ + CYr-Ys) 2 , 


« = ^(*Q,-X 6 )a + (*Q y -n) 2 , 

<j>= 1 if 6 > A 
0 otherwise 

Q* = * T X r , 

Qy = ^ T Y r , 


with 

X r = {X rl ,X r 2 ,...,X rm } T , 
Y r = {y r i,y r 2 ,...,F rm } T . 


^ = krfr' + (*, + * r )f - r c < (3.a) 

^ = ky h r c + (k g + fc r )§ - T c v c h (3.6) 

The control circuitry also includes a lead 
network. The governing equations are 


d v 


* C 1 

— - = — v 4- 
dt t‘ Vv + 


dvl 

dt 


tv,, 


(4.a) 


dv‘h 

dt 



dvl 

dt 




(4.6) 


The physical displacements of the rotor at 
the bearing and imbalance locations can be ob- 
tained using the following coordinate transfor- 
mation: 

N 

X r k = ^ ^ $ki Qxi » 

t=l 

N 

Y r* = £*kiQ„i, 

1=1 

(k = Nbi,Nb2,Ni,3,Ni m b) 

The position of the shaft at the magnetic 
bearing location is measured using horizontal 
and vertical proximity sensors. They are mod- 
elled as linear first order systems. 


dvl _ 

dt t p 


(2.a) 


Four current amplifiers (one for each coil) 
supply current to the magnetic bearing. They 
are modelled as 


dt? 

dt 


di\ 

dt 



di$ 

dt 


I , _Ll a l 

~ r a*l + T a K * V v 

(5.a) 

1 *a 1 1 

to t* fco V 

^-a * ' y-a 1 v 

(5.6) 

”1 

(5.c) 

-^•4 + 

(5 -d) 


dv* _ k v x - v v h 
dt r v 


( 2 . 6 ) 


The active control system is based upon a 
proportional-derivative (PD) control law and 


The forces supplied by the magnetic bearing 
to the rotor are assumed to be decoupled in the 
horizontal and vertical directions. The effects of 
flux saturation are included in the model based 
upon the approach described by Lewis (1993). 
The equations for these forces are 



F r ,am6 = 1*oANt( “ ]T 

(2 (D - x N b2 ) + — ) 2 

Mr 

l l - ) 

(2 (D 4- x N b3 ) + — ) 2 

Mr 

(6m) 


Fy f am 6 — fioANj* ( 


il 


( 2 (D - VNb?) + — ) 2 

Mr 


*3 


(2 {D + Vn b a) 4 ) 2 


Mr 


( 6 . 6 ) 



Figure 3 Block Diagram of Simulation Model 


DISCUSSION 

Using the rotordynamic model described in 
the preceding paragraphs, a study of the com- 
bined dynamics of a magnetic bearing sup- 
ported rotor with auxiliary bearings was per- 
formed. The parameters used in this work were 
identified from experimental evaluation of the 
test rig and are shown in Table 1. All of the 


Parameter 

Value 

Units 

K| X , K|y 

17,510 

N/m 

C.x, c. y 

2.0 

N-s/m 

K c 

87,557 

N/m 

Cc 

0 

N-s/m 

K b 

17,510 

N/m 

Cbxi Cj3y 

2.0 

N-s/m 

M b 

0.033 

kg 



kg-m 



- 


13 

“ 

k r 

23 

- 

k a 

0.45 

amp/v 

T*> 

WMESESm 

- 

r‘ 

HES23XS9I 


Tj r~ 

MirrtT-BM 

- 

T a 


- 

A 

f fTTT‘MT 

m 2 

N t 

164 

- 

D 


m 

Mo 


- 


4.2 

volt 

< 

1.9 

volt 

4' 


volt 

< 

1.9 

volt 


Table 1 Simulation Model Parameters 


responses illustrated in the figures are for this 
base parametric configuration, unless otherwise 
indicated. Figure 3 shows the variation of sys- 
tem natural frequencies with rotor speed for the 
base parametric configuration (without auxil- 
iary bearing contact). The rotor speed of 200 
rad/sec was selected as a reasonable value near 
a critical speed in order to illustrate a worst 
case behavior of the system. 

Figure 4 shows the rotor response under nor- 
mal operating conditions. The horizontal and 
vertical responses are of the same amplitude 
and frequency at steady state. The response 
consists of a transient region of about 0.5 sec- 
onds in length followed by steady synchronous 
oscillations. The initial conditions are set so 
that the rotor contacts the auxiliary bearing 
during the transient phase of the oscillations. 
The rotor vibrates downward, contacts the aux- 
iliary bearing, bounces several times, and then 
settles into a steady synchronous whirl of about 



































0.05 mm in amplitude. 

Of particular interest in this study is the use 
of the remaining coils of the AMB (when one 
fails) to provide a sideloading force to reduce 
rotor vibration. The test configuration consists 
of allowing coil 2 to fail and investigating the ef- 
fects of sideloading on the system. Coils 1 and 3 
are kept in the same electronic configuration as 
if no failure had occurred. However, the voltage 
in coil 4 is set to a value that will pull the shaft 
horizontally over to the auxiliary bearing. Two 
basic scenarios for rotor control during failure 
are considered. Figures 5. a — 5*c illustrate the 
rotor behavior under the conditions for the first 
scenario, when the rotor speed is maintained at 
a constant value. For each case, the AMB fails 
at t/= 0.5 seconds. There is some perturbation 
of the vertical motion of the rotor when the hor- 
izontal coil fails due to coupling of the horizon- 
tal and vertical motion (primarily) through the 
gyroscopic terms. However, the steady state 
vertical responses are not significantly affected. 
Figures 6. a - 6.c show the rotor behavior for 
the same failure conditions (with sideloading) 
as above for the second scenario, when the to- 
tor is decelerated at a constant rate upon AMB 
failure. Again the response amplitudes tend to 
be lowest for smaller clearance values. It is ob- 
served that the responses may be quite high as 
the rotor contacts the auxiliary bearing, even 
with sideloading from the remaining coil. How- 
ever, sideloading is shown to be quite effective 
in encouraging rotor contact that will dissipate 
the vibrational energy. The response ampli- 
tudes tend to decrease as the auxiliary bearing 
clearance is decreased. However from a practi- 
cal perspective, there are definite design trade- 
offs that must be considered. If the clearance 
is too small, the rotor may strike the auxiliary 
bearing during normal operation of the mag- 
netic bearing if the rotor is perturbed. This 
tends to unnecessarily shorten the life of the 
auxiliary bearing. If the clearance is too large, 
the rotor vibration (upon AMB failure) will be 
excessively high ana may damage the rotor or 
the magnetic bearing structure. 

The effect the sideloading voltage on the sys- 
tem dynamics is particularly interesting. The 
magnetic bearing forces are nonlinear functions 
of the dynamic clearance and the current. If one 
cod fails and the remaining coils are used to pro- 
vide sideloading, the natural approach might 
seem to issue a constant command voltage, with 
the objective of having a constant sideloading 
force at steady state. However, the inherently 
unstable behavior of the magnetic bearing with- 
out active control coupled with the vibrational 



rotor speed (rad/sec) 


Figure 3 Rotor Natural Frequencies as 
a Function of Rotor Speed 
(for the parameters of Table 1 ) 
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Figure 4 Rotor Response Without AMB Failure 
Jl = 200rad/sec, Cl = y 


effects between the auxiliary bearing and the 
rotor produces a severely unstable response. A 
remedy is to command a voltage that is pro- 
portional to the horizontal shaft displacement. 
That is, v a = v*(D + where t>“ is a con- 

stant. With an ideal amplifier, any value of v 0 
can be used with such a system. However due 
to the delay in the dynamics from the time con- 
stant of the power amplifier, unstable behavior 
may also result if t>“ is too large. Figure 7 shows 
the maximum that can be commanded with 
stable behavior as a function of power amplifier 
time constant. These results were developed by 
fitting a smooth curve through data points ob- 
tained by observing the behavior of the simula- 
tion model for various sideloading voltages and 
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Figure 5.a Rotor Response With AMB Failure 
ft = 200rad/sec, Cl — 
tj = 0.5, a = 0 



Figure 5.c Rotor Response With AMB Failure 
ft = 200rad/sec, Cl — 'y 
</ = 0.5, a = 0 
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Figure 5.b Rotor Response With AMB Failure 
ft = 200rad/sec, Cl = y 
ty = 0.5, a = 0 

amplifier time constants. It must be concluded 
that power amplifier dynamics are a critical fac- 
tor in the selection of appropriate levels of bear- 
ing sideloading. 

CONCLUSION 

A simulation model has been developed for a 
magnetic bearing supported flexible rotor sys- 
tem with auxiliary bearings including frictional 
effects. The model has been described in detail 
with experimentally obtained model parame- 
ters. The response characteristics for a vari- 
ety of system configurations were presented and 
discussed, including the effects of sideloading. 
Some guidelines were given for the selection of 
appropriate levels of sideloading. 


Figure 6 . a Rotor Response With AMB Failure 
ft = 200rad/sec, Cl = 
tf — 0.5,0; = — 100 
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ABSTRACT 

This study examines the response of a flexible ro- 
tor supported by load sharing between linear bear- 
ings and an auxiliary clearance bearing. The objec- 
tive of the work is to develop a better understand- 
ing of the dynamical behavior of a magnetic bear- 
ing supported rotor system interacting with auxil- 
iary bearings during a critical operating condition. 
Of particular interest is the effect of coupling be- 
tween the bearing/housing and shaft vibration on 
the rotordynamical responses. A simulation model 
is developed and a number of studies are performed 
for various parametric configurations. An experi- 
mental investigation is also conducted to compare 
and verify the rotordynamic behavior predicted 
by the simulation studies. A strategy for reduc- 
ing sychronous shaft vibration through appropriate 
design of coupled shaft/bearing/housing vibration 
modes is identified. The results are presented and 
discussed. 

NOMENCLATURE 

C = damping, N-sec/m 

K = stiffness, N/m 

Mb = auxiliary bearing mass, kg. 

Mh — housing mass, kg. 

N = total number of modes considered 

Nbi = node number at leftmost bearing 

Nb2 = node number at rightmost bearing 

Nb 3 = node number at auxiliary clearance bearing 

Q x = rotor modal coordinate vector in X direction 

Q v = rotor modal coordinate vector in Y direction 


t = time, s 

X r = rotor physical coordinate vector in X 
direction, m 

Y r = rotor physical coordinate vector in Y 
direction, m 

Xb = auxiliary bearing physical coordinate vector 
in X direction, m 

Yb = auxiliary bearing physical coordinate vector 
in Y direction, m 

Xh = housing physical coordinate vector in X 
direction, m 

Yb, = housing physical coordinate vector in Y 
direction 

A = radial clearance in auxiliary bearing, m 
I a = rotor polar mass inertia matrix, kg-m 2 
^ = rotor free-free modal rotation matrix 
$ = rotor free-free modal displacement matrix 

r = * T i a * 

= rotor operating speed, rad/s 
u; n = matrix of rotor free-free natural 
frequencies, rad/s 
ip = imbalance vector 
Subscripts 

bx = auxiliary bearing, x-direction 
by = auxiliary bearing, y-direction 
c = contact 

hx = housing, x-direction 
hx = housing, y-direction 
xr = rightmost bearing, x-direction 
xl = leftmost bearing, x-direction 
yr = rightmost bearing, y-direction 
yl = leftmost bearing, y-direction 
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INTRODUCTION 

In recent years, the use of active magnetic bear- 
ings for supporting turbomachinery has been an 
area of interest for both academic researchers and 
turbomachinery users in industry. Magnetic bear- 
ings provide the potential for significant improve- 
ments over other types of rotor support, including 
elimination of wear and bearing friction-related en- 
ergy losses as well as a means of actively suppress- 
ing rotor vibration. However, their use has been 
significantly limited due to a number of technical 
problems. A particular area of concern is the aux- 
iliary bearing, which protects the soft iron core of 
the magnetic bearing and provides rotor support in 
case of overload or failure of the magnetic bearing. 

Typically, the auxiliary bearings have relatively 
small clearances so that magnetic bearing ro- 
tor/stator contact does not occur during bearing 
failure or power loss. Due to these small clear- 
ances, contact between the rotor system and the 
auxiliary bearings can occur during standard oper- 
ation of the magnetic bearings. When this occurs, 
load sharing between the magnetic bearings and the 
auxiliary bearings results, and the rotor system in- 
teracts with its auxiliary bearings. The dynamics 
of such an occurrence must be understood in order 
to properly design the auxiliary bearing system for 
a magnetic bearing supported rotor. 

A number of different bearing types have been 
suggested as auxiliary bearings. These include 
bushings, rolling element bearings, and various 
types of journal bearings. The most commonly 
considered are rolling element bearings. The ma- 
jor disadvantage associated with using rolling ele- 
ment bearings (or bushings) is the requirement of 
a clearance between the rotor and the inner race of 
the bearing, without which many of the advantages 
associated with using magnetic bearings would be 
reduced or eliminated. This clearance introduces 
a nonlinear dynamical feature which may signifi- 
cantly impact the behavior of the rotor. 

There are quite a number of studies in the lit- 
erature concerned with nonlinear rotordynamics. 
Ehrich (1965 and 1967) studied the rotor/stator 
interaction with a clearance and found zones of 
bistable synchronous behavior. Black (1968) also 
found these bistable interaction zones. Ehrich ex- 
tended this research by predicting subharmonic 
behavior (1966 and 1988), superharmonic behav- 
ior (1992), and chaotic responses (19911. Bently 
(1974), Muszynska (1984), and Childs (1979 and 
1982) also observed and studied subharmonic re- 
sponses due to clearance effects. 

There is relatively little work available in the 
open literature that is specifically concerned with 
the dynamics between the rotor system and the 
auxiliary bearings. The research that has been 


performed to date is primarily concerned with the 
transient behavior of rotor drop on the auxiliary 
bearings due to power failure or inoperative mag- 
netic bearings. Two papers with this focus are 
Gelin et al., (1990) and Ishii and Kirk (1991). 

The present work is concerned with developing 
an understanding of the dynamic behavior of a 
rotor system supported by load sharing between 
the linear bearings and an auxiliary bearing with 
clearance. Of particular interest is the influence of 
coupled shaft /bearing/housing vibration modes on 
the rotordynamical behavior of such a system. 

SIMULATION MODEL 

Figure 1 shows the simulation model used for 
this investigation. It consists of a flexible rotor 
supported at both ends by magnetic bearings. A 
rigid disk with adjustable imbalance is placed at the 
midpoint of the bearing span. As a simplification, 
the magnetic bearings are modeled as spring and 
dashpot systems and interaction with a single 
auxiliary bearing is considered. The auxiliary 
bearing was modeled as an antifriction bearing with 
a clearance and a mass attached to a housing. The 
housing in turn has a mass, stiffness, and damping. 
Figure 2 shows the auxiliary bearing model used in 
the simulation. 


Rigid Disk 



Magnetic Bearing Magnetic Bearing 


Figure 1 Schematic Diagram of Simulation Model 

The rotor is modelled using the free-free bend- 
ing mode shapes and natural frequencies obtained 
through finite element analysis. The finite element 
code uses 19 stations and the first four modes (two 
rigid body and two flexible modes) are included in 
the simulation model. 

Using these simplifications, the equations of mo- 
tion for the system can be written as follows: 

Q r +OrQ y + u;’Q r + $ T (F x , /in + F r , n/ ) 


2 



= * T F x , imb , (l.a) 

Qv-ftrQx + <4Q„ + * T (F yMn + F v< „,) 


— F 

(1.6) 

M b X b +^(X 6 - X b ) + K bx (X b - X h ) 


— F t^nh 

(l.c) 

M b Y b +Of-(Y b -Y b ) + K by (Y b - Y h ) 


= Fy t nl > 

(Yd) 

M h X h +( ^ + (K bx + K hx )X h 

= C bx X 6 + K bx X bl 

(Ye) 

M H Y h +ty + + (A- k „ + K h ,)Y h 


= 

(Yf) 


where 


f x ,k„ = * n »(k x1 q x + -fq x ) 

+ *N„(KttQt + %"Q*) 

F Vt/in = + ^fQ*) 

+ ^ M (ifyrQv+% 1 Qv) 

F x , imt = xl>n 2 cosSlt, 

Fyjmb = i>n?sinnt, 

F Ml = ^-A)( §Q y^ ) 

Fy, n / = (j>F conii^ - A)( #Qy 5 ~ — ) 

^y^-x^ + ^Q.-n) 2 , 

<f> = 1 if 6 > A 

0 otherwise 
Qx = * - 1 X r , 

Qy = ^ _1 Y r , 


with 


X r — {X r i, X r 2 J . * • , Xrm} > 
Y r = { Y rU Y r2 ,...,Y rm ) T • 

(m = total number of nodes) 


The physical displacements of the rotor at the 
auxiliary bearing or magnetic bearing locations 
can be obtained using the following coordinate 
transformation: 


N 

Xrk = ) ^ $ki Qxii 

i= J (k = N bl ,N n ,N b 3 ) 

Yrk = ^kiQyi, 

1 = 1 


Rotor 



Figure 2 Auxiliary Bearing Model 


EXPERIMENTAL MODEL 

Experimental tests were performed in order to 
validate the behavior predicted by the simulation 
model and to gain some insight into the dynamical 
responses that are to be expected. A schematic 
diagram of the rotor test rig that was used in the 
experimental work is shown in Figure 3. 

The rotor used in this study has two basic 
components: a flexible shaft and an auxiliary 

clearance bearing. The shaft is made of steel 
and is 0.374 inches in diameter and 18.0 inches 
in length. It is supported at 1.0 inch from the 
right end by ball bearings suspended in a frame 
by four springs and at 1.0 inch from the left end by 
a bushing with a tight clearance. These supports 
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Figure 3 Experimental Model 


represent the magnetic bearings. The stiffness of 
the left support is 17511 N/m (101.5 lb/in) for 
both the horizontal and vertical directions. It 
is used to somewhat isolate the rotor from the 
effects of the flexible coupling which attaches the 
rotor to the motor and to enforce low amplitude 
vibration at this location to protect the motor. The 
stiffness of the right support, for both horizontal 
and vertical directions, is 2539 N/m (14.5 lb /in). 
This lower stiffness allows for significant vibration 
of the rotor in the speed range of the motor. A 
rigid disk with an adjustable imbalance is placed 
at the midpoint of the bearing span. The auxiliary 
clearance bearing/housing consists of a bushing 
suspended in a frame by four springs (Figure 4). 
It is situated at the right end of the rotor. The 
clearance is adjustable by changing the bushing. 
The auxiliary bearing/housing stiffnesses are varied 
by interchanging the springs. 

The rotor is driven by an adjustable speed mo- 
tor with feedback speed controller. Shaft vibration 
is measured using eddy current proximity displace- 
ment sensors fixed so as to measure displacement 
in the vertical and horizontal directions. The dis- 
placement signals were recorded and analyzed with 
a signal analyzer. 


Frame 



DISCUSSION 

A series of parallel studies were conducted using 
the simulation and experimental models discussed 
in the preceding sections. The following discussion 
is aimed at comparing the predictions of the simula- 
tion model with experimentally observed responses, 
with the objective of obtaining insight into the be- 
havior of flexible rotor system dynamics due to load 
sharing between the magnetic bearings and an aux- 
iliary bearing. 

Using the data from measurements and calcula- 
tions, the stiffness, damping, and mass character- 
istics of the experimental rotor system were identi- 
fied. The imbalance configuration consists of 3.0 x 
10“ 6 kg-m located at the left end of the rotor shaft 
(from the flexible coupling to the electric motor) 
and a variable imbalance located on the rigid disk. 
Linear, hysteretic, and coulomb friction damping 
models were considered. The support damping ap- 
pears to result principally from hysteresis. If syn- 
chronous vibration is assumed, this results in a 
damping coefficient scaled by the rotor speed. The 
numerical values are given in Table 1. 

The governing parameters for an auxiliary bear- 
ing are stiffness, mass, damping, clearance, and ax- 
ial location. The axial location generally is fixed 
by other considerations, such as space requirements 
and the need to be close to the magnetic bearing 
to better protect it. This study focusses on the 
remaining design parameters. For the system stud- 
ied, significant nonsynchronous vibration occured 
only for the 2 u component and for a limited ro- 
tor speed range. Figure 5 illustrates the observed 
behavior for various imbalance values. The ampli- 
tudes are quite small and, interestingly, virtually 
independent of imbalance. The remainder of this 
discussion will focus on the synchronous response 
amplitudes. 
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Parameter 

Value 

Units 

K b 

471 

N/m 

K h 

1.8xl0 7 

N/m 

K c 

87,557 

N/m 

K| X , Kly 

17,510 

N/m 1 

K„, K ry 

2,539 

N/m 

C| X ? Qy 

1,000 

N/m 

CfX? Cry 

500 1 

N/m 

C|y, C|y 

1,000 

N/m 

Cry * C ry 

500 

N/m 

Cbx? Cby 

2,000 

N/m 

Chx> Chy 

0 

N/m 1 

M b 

0.033 

kg 

M h 

0.296 

kg 

01 

3.0 x 10~ 6 

kg-m 

V’lO 

2.5 x 10“ s 

kg-m 

A 

5.0xl0 _i 

m 


Table 1 Simulation Model Parameters 



Figure 5 Experimental Results for 
2u> Component 

The simulation responses discussed below were 
determined from a single-term harmonic balance 
analysis that was numerically implemented and ver- 
ified at selected points using direct numerical inte- 
gration of the governing equations of motion. The 
harmonic balance procedure is described by Nayfeh 
and Mook (1979). Figures 6.a - 6.e show the re- 
sponse amplitudes predicted from the simulation 
study for the parameters of Table 1. Figures 7. a - 
7.c show some corresponding experimental results. 
The rotor speed range extends to above the second 
critical speed, with the first critical speed primar- 


ily a rotational mode and the second critical speed 
primarily a translational mode. 

Figures 6. a presents simulation results for vary- 
ing auxiliary bearing stiffness and Figure 7. a shows 
the corresponding experimental values. There is 
relatively good agreement between the actual and 
predicted responses except for the peak values 
which are very sensitive to damping. At lower ro- 
tor speeds, the stiffness effect from the auxiliary 
bearing tends to dominate and serves to increase 
the effective lower critical speed. At higher ro- 
tor speeds, the inertial effects of the bearing tend 
to dominate and the added mass serves to lower 
the effective second critical speed. This effect is 
balanced somewhat for higher stiffnesses (as is ex- 
pected), as can be seen for the case with Kb =2,625 
N /m. The damping added to the system from the 
auxiliary bearing serves to attenuate the peak re- 
sponse amplitudes for both critical speeds. The re- 
sponse amplitudes between the two response peaks 
are affected little by the auxiliary bearing, with the 
rotor and bearing/housing vibrations tending to 
decouple for response amplitudes below the clear- 
ance. Above the second critical speed, the pres- 
ence of the auxiliary bearing serves to reduce the 
response amplitudes for the remainder of the plot. 
For this region, there is significant coupling be- 
tween the rotor/bearing/housing vibration and the 
response amplitudes are below the clearance for 
speeds greater than about 550 rad/sec. The ro- 
tor speed ranges for which coupling of the ro- 
tor/bearing/housing vibration will occur depend 
upon the phase of the combined system response, 
as described by Black (1968). For the purposes of 
this study, Kb=471 N/m is chosen as a reasonable 
value and variations of other parameters are per- 
formed with this auxiliary bearing stiffness. 



w (rad/s) 


Figure 6. a Simulation Results for Varying 
Auxiliary Bearing Stiffness 
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Figure 6.b Simulation Results for Varying 
Auxiliary Bearing Mass 


Figure 6.d Simulation Results for Varying 
Imbalance 



- Figure 6.c Simulation Results for Varying 
^ Auxiliary Bearing Damping 

Figures 6.b shows the predicted rotoT response 
for the parameters of Table 1 with the bearing mass 
varied. These results show that the more massive 
bearing/housing configurations can result in much 
lower amplitude rotor responses. In fact, for suffi- 
ciently high mass values, the response amplitude is 
dramatically reduced for the rotor speed range from 
about 600 rad/sec to about 780 rad/sec where it 
merges with other amplitude response curves. Fig- 
ure 6-c shows the effect of auxiliary bearing damp- 
ing on the rotor response. While the response am- 
plitudes remain relatively unchanged for the ma- 
jority of rotor speeds, higher damping does serve 
to significantly reduce the rotor response ampli- 
tudes through the critical speeds. Damping also 
effects the phase characteristics of the combined 
rotor /bearing/housing system and changes the ro- 



Figure 6.e Simulation Results for Varying 
Auxiliary Bearing Clearance 

tor speed ranges for which coupled responses occur. 
Figures 6.d and 7.b show the rotor response for the 
parameters of Table 1 with the imbalance varied. 
As expected, the response amplitudes increase with 
increasing imbalance. Figures 6.e and 7.c show the 
rotor response amplitudes with the auxiliary bear- 
ing clearance varied. As the clearance increases, the 
dynamics of the combined rotor/bearing/housing 
vibration tend to decouple over a wider speed range 
between the two critical speeds. Again, such be- 
havior is expected from the analysis work for a ro- 
tor operating in an annular clearance performed by 
Black (1968). 

Close examination of the above results reveals a 
very interesting trend. Turbomachine are designed 
to operate at rotor speeds between the various sys- 
tem critical speeds so as to minimize synchronous 
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Figure 7.a Experimental Results for Varying 
Auxiliary Bearing Stiffness 



w (rad/s) 


-Figure 7.b Experimental Results for Varying 
Imbalance 

response amplitudes. For the present system, the 
lowest amplitude responses (for operation above 
the second critical speed) occur for the auxiliary 
bearing configurations where there is the largest de- 
crease of the second effective critical speed, which 
is due to a coupling of rotor/stator /housing vibra- 
tion. While not really suprising, this result indi- 
cates an interesting strategy for the design of aux- 
iliary support using clearance bearings. Selection 
of combined shaft /auxiliary bearing/housing vibra- 
tion characteristics such that the effective critical 
speeds nearest the operating speed (directly above 
and below) shift away from it should result in lower 
amplitude rotor vibration and provide the best op- 
portunity to protect the magnetic bearing. How- 
ever, the phase characteristics of the combined sys- 
tem must be such that vibration coupling occurs 



Figure 7.c Experimental Results for Varying 
Auxiliary Bearing Clearance 

and care must be taken to ensure that this is the 
case if the maximum benefit is to be derived from 
the auxiliary bearing. 

CONCLUSIONS 

A study of the dynamical behavior of a flexible 
rotor supported by linear bearings (representing a 
set of magnetic bearings) and an auxiliary bearing 
with clearance has been presented. Parallel sim- 
ulation and experimental studies have been per- 
formed for a variety of parametric configurations. 
The influence of bearing mass, stiffness, damping, 
clearance, and imbalance were examined. It was 
observed that the dynamic behavior of a rotor in- 
teracting with an auxiliary bearing depends very 
strongly on the structural parameters of the auxil- 
iary bearing and associated housing. Appropriate 
selection of these parameters is critical if acceptable 
vibration characteristics are to be obtained for such 
systems. In general, one must carefully consider the 
influence of housing/bearing dynamics and how to 
best take advantage of favorable coupled modes of 
vibration in the development of auxiliary bearing 
designs. 
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